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Abstract 



We present a study of phase separation from supersaturated 3 He- 4 He liq- 
uid mixtures at low temperatures addressing both the degree of critical super- 
saturation Ax cr and the thermal-to-quantum crossover temperature T* for the 
nucleation process. Two different nucleation seeds are investigated, namely 
3 He droplets and 4 He vortex lines with cores filled with 3 He. We have found 
that the experimental T* is reproduced when we consider that nucleation pro- 
ceeds from 3 He droplets, whereas Ax cr is reproduced when we consider 4 He 
vortex lines filled with 3 He. However, neither nucleation configuration is able 
to simultaneously reproduce the current experimental information on Ax cr 
and T*. 
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I. INTRODUCTION 



Supersaturated 3 He- 4 He liquid mixtures have been recognized as an ideal system to study 
nucleation phenomena. On the one hand, these mixtures can be made extremely pure, and 
on the other hand, 4 He covers the walls of the experimental cells, preventing undesired 
heterogeneous nucleation that is difficult to avoid in other substances. Moreover, the fact 
that helium remains liquid up to zero temperature (T) has raised expectations that the 
transition from thermal to quantum nucleation could be observed as T decreases. 

The first experiments on metastable helium mixtures were carried out in the sixties^!, 
but one had to wait more than twenty years for a systematic study of this phenomenon! - !. 
Meanwhile, the kinetics of nucleation in dilute helium mixtures was theoretically addressed 
by Lifshitz et al0, much along the line of the classical work of Lifshitz and Kagan0. 

At very low temperature and at zero pressure, the concentration of 3 He in the mixture 
at saturation is! x s ~ 6.6%. However, supersaturated 3 He- 4 He mixtures can be found in 
a metastable state for concentrations x above the saturation value. A first estimate of 
the degree of supersaturation Ax cr = x — x s resulted from extrapolation of the measured 
3 He chemical potential excess A/X3 along the demixing line. This extrapolation yielded 
dAfi 3 /dx > up to values of the 3 He concentration x > 16 %, thus givingEl Ax cr ~ 10%. 
Large values of critical supersaturation have also been obtained in microscopidl! and density 
functional! calculations. 

An intriguing observation is the small degree of critical supersaturation attained in recent 
experiments on phase separation in supersaturated helium mixtures. It was found to be of 
the order of 1 % in the experiments carried out by the Ukrainian groupSi'Si, and below 0.5 
% in those of the Japanese group!'!!. As a matter of fact, the measurements in the sixties!! 
were also yielding small Ax cr . 

A possible explanation of the discrepancy between theory and experiment is that phase 
separation at small supersaturations may be triggered by heterogeneous nucleation on vortex 
lines. Indeed, at the low temperatures involved in the experiments 4 He is still superfluid 
due to the limited solubility of 3 He into 4 He. This mechanism was proposed in Ref. [15|, 
and it has been further elaborated^!!. It turns out that the presence of vortex lines in the 
mixture decreases the degree of critical supersaturation from a value Ax cr ~ 10% to ~ 1 
%. At Dresent, both experimental groups seem to have adhere to this explanation of their 



result st5rE3. 

Besides critical supersaturation, another experimental quantity of most interest is the 
temperature at which the nucleation process changes from being a thermal one to a quantal 
one. This crossover temperature T* has been determined by measuring the temperature 
below which the degree of critical supersaturation becomes almost T independents; it is of 
the order of 10 mK. Obviously, the theoretical description of the process cannot be considered 
as achieved until Ax cr and T* are satisfactorily obtained within the same framework. 

The present aim is to complete our previous works with the calculation of the crossover 
temperature using the same basic ingredients that we have used to obtain the critical super- 
saturation. To this end, we have resorted to a dynamical method, namely the functional- 
integral approach (FIA), that proved^ to be well suited to quantitatively describe quantum 
cavitation in liquid 4 He. 
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II. FORMALISM 



The segregation of 3 He from a dilute 3 He- 4 He liquid mixture at low T is a first order 
phase transition. This means that in order to nucleate the pure 3 He phase, the system has 
to overcome an energy barrier AQ. Thus, the mixture may be at 3 He concentrations above 
the one corresponding to the two-phase equilibrium value at the given pressure (P) and 
temperatureil. 

At high enough T, segregation is thermally activated and the energy barrier is overcome 
by the energy provided to the system by a thermal bath. In this case, the nucleation rate, 
i.e., the number of critical nucleation seeds formed per unit volume and time is given by 



Jt = Jot exp 



k B T 



(1) 



where the prefactor Jot depends on the dynamics of the nucleation process, ks is the Boltz- 
mann's constant and AQ max is the height of the energy barrier. At low T, segregation by 
thermal activation is no longer possible. However, it can proceed by quantum tunneling: 
the metastable liquid mixture 'tunnels' through the energy barrier. In the limit of zero tem- 
perature, the transition is purely quantal, but below T* thermally assisted tunneling also 
occurs. This is one of the many thermally assisted quantum tunneling processes occuring in 
physics (see for example Ref. |23| and refs. therein). For T < T* the tunneling rate is 

J Q = J 0Q exp(-S Q ), (2) 

where exp(-S^) is the tunneling probability, and the prefactor J q is of the order of the 
number of nucleation sites per unit volume times an attempting frequency. 

To obtain the tunneling probability one formally starts from the statistical average of 
the transition probability over a time t = tf — tf 



P 



Y,ij \{^f\exp [-± £ dtii] l^expi-Ei/kBT) 



J2iexp(-Ei/k B T) 



(3) 



where 7i is the Hamiltonian of the system, E\ are the energy states, and and are 
the wave functions of the initial and final states. This expression can be written in a 
more workable form by using the path integral formulation of Quantum mechanics, and its 
connection with statistical mechanics. According to FeynmanEa (see also Refs. |25|-f27|) , P 
can be written as the functional integral 



D[q{r)] exp 



- f drC[q(r)]\ , (4) 



where £[g(r)] is the imaginary-time (r = it) classical Lagrangian of the system and D[q(r)] 
denotes integration over all periodic trajectories q(r) with period t p = h/ksT in the potential 
well that results from inverting the energy barrier. The integral in the exponent of Eq. (|J) 
is the imaginary-time action S(T) = § drC evaluated over the period r p . In the semiclassical 
limit S(T) » h, the trajectory that contributes the most from all possible periodic orbits 
is the one which minimizes the action. It leads to Eq. (|2|) with = S m i n (T) /h, where S m i n 
is the minimum action. 
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The practical usefulness of Eq. (|J) for the problem at hand is that physical insight 
allows one to guess £ as a functional of densities and collective velocities instead of dealing 
with the impracticable Eq. (|3]). To implement this scheme, it is quite clear that a sound 
approximation is needed for the imaginary-time action. This means a realistic energy barrier 
and a simple yet reliable choice of the integration path q(r). 

Let us now work out in some detail the case in which only one collective coordinate 5 is 
considered. It is then a simple task to minimize the imaginary-time effective action: 



S(T) = idrC[8{r)} = dr 

J J-Tp/2 



M(5)5 2 + AQ(5) 



(5) 



As indicated, the effect of continuing the action to imaginary time is to invert the 'potential', 
i.e., AQ — > —AQ in the equation of motion, and the identification of /c#T with h/r p ; the path 
8(t) defined in imaginary time r has to fulfill the periodic boundary condition 5(—t p /2) = 
S(t p /2). This is illustrated in the schematic Figure |l| We have supposed that the collective 
mass M depends on 5, which will be the practical case. Imposing the extremum condition 
on the action yields the following equation of motion for S(t): 



IdM - 2 dAQ 



(6) 



Multiplying Eq. (^) by 6 we have 

d 

d^ 



-M{5)6 2 - AQ(6) 



. 



Thus 



with AQ r 



^M(5)5 2 - AQ(8) = constant = -E 



(7) 



(8) 



_ > E > 0. 

Eq. (13) has the trivial solution S = 5q corresponding to the minimum of — AQ: the system 
is 'at rest' at the bottom of the inverted barrier potential well. In this case, E = AQ max , 
and the integration of Eq. (H) yields 



1 



Aa 



r p /2 
-t p /2 



dr 



AQ r 



(9) 



Thus, the trivial solution yields the exponent for classical thermal activation, see Eq. ([[]). 
It means that within FIA, the transition between thermal and quantal regime is smooth. 
For E < AQ max , one has to seek periodic solutions S(r) whose turning points 5i and 5 2 are 
such that AQ(Sx) = AQ(5 2 ) = E (see Fig. p. Integrating Eq. @ we get the period r p : 



, . rHE) M(5) 



(10) 



Using 



4 



dr 

the action Eq. @ becomes 

2 r^(E) 



s °' £ >— £,^ |2AfiW - £ H 2[An ( l)U] - (12) 

At T = 0, t p = oo and E = 0. In this case, the solution to Eq. (|8|) is the usual instantonll, 
and S®(E = 0) coincides with the WKB approximation at zero energj 



s wkb = d5 J2M(5)AQ(5) . (13) 

hJ5i(E=0) V 

The crossover temperature is obtained by equating Eqs. ([I]) and (0) we get: 

S^(E = AQ max ) = . (14) 

Taking a trajectory 5 corresponding to E w Af2 ma:r and using Eq. (|T3), we can write 

2 fl52(i?RiA!J maI ) 



_ 2 fOj^ramifflMj M(S) 

v ma.; h ™*J Si(E ^ AUmax) \2[An(5)-E] 



(15) 



The comparison to Eq. fllOD yields ^/(/^T*) = r p (E pa AQ max ). An analytical expression 
for T* is obtained by expanding Eq. (|i~0|) around the maximum of AQ located at Sq. It 
reads 



h 


1 


d 2 Att 




2n\ 


M(S ) 


dS 2 


So 



(16) 

Eq. ( |T6| ) shows that the value of T* is determined by small variations around 5q. This is a 
well-known result that was suggested long ago by Goldanskiiil, the validity of which goes 
beyond the simple model we have used to derive it. Generally speaking, the crossover temper- 
ature is determined from the frequency of the small amplitude oscillations around the min- 
imum of the inverted barrier potential well. When more realistic methods are employed to 
generate the critical nucleation configurations, like for example density functional theorv^S, 
the problem becomes an infinite-dimensional one, whose solution is quite a formidable task. 
Crossover temperatures for thermally assisted quantum cavitation in liquid helium have 
been determined only recently0'@. 

Finally, we mention that S WKB can be used to estimate T* through the expressionSS 

1 rp* _ AVL max 

KB J- WKB — gWKB ' v '/ 
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III. ENERGY BARRIERS AND IMAGINARY-TIME ACTIONS 



As indicated at the introduction, we have considered two different nucleation seeds, de- 
pending on whether there is some vorticity in the experimental sample or not. The first 
one is a pure 3 He droplet embedded into the 3 He- 4 He mixture, and the second one is a 4 He 
vortex line filled with pure 3 He. In the following we work out in detail the energy barriers 
and imaginary time actions making use of the capillarity approximation in the case of a 
droplet, and of the hollow core approximation in the case of a vortex line. These are sound 
approximations provided nucleation occurs near the two phase equilibrium curve. Indeed, it 
is well known that near this curve, the critical nucleation configuration, i.e. that correspond- 
ing to the maximum of the energy barrier, is large. This makes negligible 'finite size effects' 
such as interface diffuseness and curvature corrections. This has been checked by comparing 
the capillarity model results with density functional calculations, which selfconsistently in- 
corporate finite size effects§ii. We recall that experimentally, the process occurs near the 
saturation line; otherwise, Ax cr would be large. This is at variance with the process of cav- 
itation in liquid helium, which takes place near the spinodal lineEH 3 !, where the cavitation 
seeds are small, pachydermic bubbles and the capillarity approximation fails§ffl. In the 
case of nucleation triggered by 3 He droplets, we will also present the results obtained using 
a density functional FIA as described in Ref. |32|. 



A. 3 He droplets 

Within the capillarity approximation, the energy barrier of a 3 He droplet of radius R 
immersed in the mixture is written as 

Atl(R) = A-kR 2 <t - —R 3 p 30 A/13 , (18) 

where a is the surface tension of the 3 He- 4 He interface, P30 is the density of pure 3 He inside 
the droplet, and A/X3 is the difference between the chemical potential of 3 He in the metastable 
mixture and in pure 3 He. This difference is negative when the mixture is stable and positive 
otherwise. Upon minimization with respect to R, the maximum of the energy barrier and 
the radius of the critical droplet R c are given by 

AO 47T 2 167T (X 3 



3 c 3 (p3oA/! 3 ) 2 
R c = 2a/{ Pm A^) . (19) 

Values of R c and AQ max can be obtained using the experimental valuesBflfi of the 
physical quantities entering their definition. In particular, the experimental data@ can 
be employed to write the above equations as a function of Ax using the approximation 
A/i3 ~ (dAfis/dx)Ax along the demixing line. For instance, at zero pressure A/13 ~ 2.49 Aa; 
K. In view of the weak T dependence of these quantities in the range of temperatures of 
experimental interest (below 150 mK), it is legitimate to consider their values at zero K. 
If nucleation proceeds thermally, the degree of critical supersaturation at fixed P and T is 
obtained by solving the equation 1 = (V ■ t) exp JT- It determines Ax cr from AQ max (Eq. 
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AQ max kg 

T\n[(V-t) exp J 0T }. (20) 

We have used as prefactor an attempting frequency per unit volume J t = Vq/Vq ~ 
{k B T /h) / (AhRq/?>). As pointed out, the prefactor depends on the physics of the nucle- 
ation process, but the results are not very sensitive to the precise value of the product 
iy ■ t)ex P JoT (see the discussion after Eq. (^)). Yet, using the estimates made at the end 
of Ref. |32] we have checked that in the present case, where phase separation proceeds by the 
diffusion of 3 He atoms, the diffusion time to make a critical 3 He droplet of Rq ~ 10A is, in 
the worst case, of the same order of magnitude than the 'thermal period' tt = 2irh/kBT. 
As thoroughly discussed by Burmistrov et alH, care should be exerted when describing later 
stages in the dynamical evolution of the critical droplet after it reaches a mesoscopic size. 
Taking for the experimental volume and time (V ■ t) exp those of Ref. ^5| and for R ~ 10A, 
we obtain Ax cr ~ 13% at T = 100 mK, and ~ 19% at T = 50 mK. This result is an order 
of magnitude larger than the experimental data. This is a drawback of using 3 He droplets 



as nucleations seeds detected inmediately after the first systematic experiments 

Let us now calculate the crossover temperature. In the following, we write R = R c + 5, 
where here 5 measures the displacement of the surface of the drop from the value corre- 
sponding to the critical configuration. In this way, the dynamics of the nucleation process 
is described by the time evolution of 5. The energy barrier then reads 

Air 

AQ(5) = MR C + 5fa - —{R c + <5) 3 p 3 o A/z 3 • (21) 

The next step is to express the imaginary-time Lagrangian in terms of the displacement 5 
and its time derivative 5. Since the 'potential part' is just — AQ(5), we have only to derive 
the kinetic energy term. This last is given by 

E kin = ~Y J d rp3(r,t)u^(f,t) + ^ J drp 4 (r,t)u%(r,t) 

= l -M{5) 5 2 , (22) 

where p q (r,t) and u q (r,t) are respectively, the particle density and the collective velocity 
of isotope q. The collective mass M(5) is determined as follows. For a given configuration 
characterized by a 5 value, the densities of the system are written as: 

pi(r) = p 45 (r) = p 4m [l - Q(R C + S-r)] = p ic {r - 5) 

p 3 (r) = p 35 (r) = p 3m [l - Q(R C + 5 - r)] + p 3 o©(^ c + 6 - r) = p 3c (r - 5) , (23) 

where p qm with q = 3,4 is the particle density of each isotope in the metastable phase 
(x = P3m/(p3m + Pim)), and the subscript c refers to the critical configuration. This is 
schematically illustrated in Fig. [| The dynamics comes into Eqs. (p3|) through the time 
dependence of 5, i.e., p q (r,t) = p qc {r — 5{t)). Thus, 

p q (r,t) = -5-^p qc (r) . (24) 

Assuming spherical symmetry during the process of growing, the densities and collective 
velocities only depend on the modulus of r, and it is possible to write u q (r,t) as a function 
of 5 and 5 formally integrating the continuity equation^ 
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^ + V( Pq u q ) = 0. (25) 



This yields 



u q (r,t) = — - / s p q {s,t)ds . (26) 

r 2 p q (r,t) Jo 



Using Eqs. (^f) and (|2~6"D the mass parameter becomes 

(P30 - P3r. 



M(5) = 4tt 



m 4 p4m + m 3 - 



(Rc + 



P3m 

= AnM(R c + 5) 3 . (27) 

The same expression has been derived in Refs. |TB1 , |3"3"| using a different method. 

By inserting the expressions for Afi(5) and M(5) in Eq. ([16]), we straightforwardly 
obtain the crossover temperature 



Furthermore, Eq. fllTf) yields 




T wkb = \I^^T*^ 1.03 1*. (29) 



^From the Ax cr values quoted below Eq. ([20]) , we estimate T* ~ 16 mK for Ax cr = 13 %, 
and T* ~ 28 mK for Ax cr = 19 %. These crossover temperatures are compatible with those 
experimentally foundi. We want to stress that these estimates are quantitative provided 
Ax is small, i.e.^far enough from the spinodal line in the P — x plane, which is determined 
by the conditioro (d[i 3 /dx)p j T = 0. Indeed, as the nucleation barrier approaches to zero at 
the spinodal line, so does T*, and the law T* ~ (Ax) 3 / 2 following from Eqs. fll9h a nd (p8|) 
eventually breaks down. A similar breakdown occurs in the cavitation processt^KLl. Fig. |3| 
displays T*(Ax) obtained from equation (p8|) at P = and 1 bar (dashed lines). The values 



of P3 m and p4 m entering Eq. (P7|) when x is much larger than x s have been calculated by 



means of the density functional of Ref. Had we used the values at the saturation line, 
the resulting T* would have been a factor of two smaller for Ax larger than ~ 5%. 

We have also obtained T*(Ax) using density functional theory to describe the critical 
nucleation clusters, and the FIA for the nucleation dynamicsil. The results are displayed 
in Fig. [3] (solid lines). For completeness, and to clarify some aspects underlying our way of 
obtaining the collective mass in the capillarity approximation, we present in the Appendix 
the equations of motion used to determine T*. 

It can be concluded from Fig. |3| that the T*(Ax) function obtained using a density 
functional that accurately describes the thermodynamical properties of the liquid mixturee3 
has a maximum of ~ 25 mK at Ax ~ 18 %. Consequently, at P = above T ~ 25 mK the 
nucleation of the 3 He-rich phase always proceeds thermally. At variance with the capillarity 
method, the density functional approach is able to realistically describe critical nucleation 
configurations from the saturation to the spinodal point. It is worth to notice that the 
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crossover temperatures obtained using the density functional plus FIA tend to zero at the 
saturation and spinodal points, whereas those obtained from the capillarity approximation 
vanish only at saturation. The calculated spinodal Ax value at P ~ is ~ 25% (see Ref. P 



and also Fig. |3|). The crudeness of the capillarity approximation when applied too far from 
the coexistence line can be assessed comparing the density profiles of the critical nucleation 
configuration of the density functional method versus the capillarity approximation. Fig. |] 
shows them at P = for Ax cr = 16 %. 

Figure ^| simply tells us that if the mixture is brought to a supersaturation degree Ax, 
nucleation will proceed by quantum tunneling below T*(Ax). To determine which of the 
T*(Ax) would correspond to a given experimental situation, one has to determine the critical 
degree of supersaturation Ax cr , which corresponds to the largest one the mixture can sustain 
before 3 He droplets nucleate at an appreciable rate. This is obtained by solving the equation 

1 = {V ■ t) exp J (30) 

taking either J = Jt for T > T* or J = Jq for T < T*, and using as a prefactor J q the 
estimate k B T*/(hV ), where V = 4nR%/3 with R = 10A. Within the more accurate density 

2 3 

functional FIA and using as experimental (V ■ t) exp two extreme values, namely 10 14 A sec 
and 10 4 A 3 sec, at P = one gets Ax cr ~ 18% and T* ~ 25 mK, and Ax cr ~ 16% and T* ~ 
23 mK respectively, and very similar values at P = 1 bar. The maximum of the thermal 
nucleation barriers AQ max are ~ IK and ~ 0.5 K, respectively. 



B. 4 He vortex lines filled with 3 He 



The possibility of considering vortex lines as seeds of a kind of heterogeneous nucleation 
in the mixture stems from the fact that the experimental sample, if treated in conventional 
fashion, is likely permeated ab initio by quantized vortices stabilized by surface pinning^. 
These vortex lines, whose core is filled with 3 He, are stable for x < x«, and metastable 



otherwisellSO. Within the hollow core model (HCM), is it very simple to obtain their 
structure. The energy per unit length of vortex line is 



tt(R) = 2nRa 



h 2 

nR 2 pw A/j. 3 + Ten 2 — p 4ri 
m 4 



In 



Roc 



(31) 



In this equation, R is the radius of the vortex core, R^ is a large enough radius at which 
velocity vanishes, n — 1,2,... is the quantum circulation number, and the remaining vari- 
ables have the same meaning as before. Minimization with respect to R yields at most two 
extrema corresponding to the values 



R — 2i?o 



/'c 



ld= 4/1 — 



A/i 3 



where 



R = 



He = 



n 2 Ti 2 p im 
2o"m 4 
<r 2 m 4 
2n 2 h 2 p 30 p 4m 



(32) 



(33) 
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Ro = limA/ig^o-R is the radius of the vortex core at the demixing line, i.e., A/i 3 = 0. If the 
vortex is in the stable region of the P — x plane, A/i 3 is negative, only the minus sign in Eq. 
(|32|) makes sense and the corresponding R defines the radius of the stable vortex. If A// 3 > 0, 
Eq. (|32|) yields two core radii. The smallest one P< corresponds to a minimum of Q(R), 
and the largest one P> to a maximum. This is, of course, the situation for given (P, T, Ax) 
values. Consequently, ab initio stable vortex lines become metastable as the mixture is driven 
into the supersaturated region. These metastable vortex lines can trigger phase separation 
by thermal fluctuations if an energy per unit length equal to AQ max = [f2(P>) — fi(P<)] 
is subministrated to them. As before, calling R c the radius P> of the critical configuration 
and defining R = R c + 5, the energy barrier per unit length reads 



AQ{6) = 2n(R c + 5)o - tt(P c + 5) 2 p 30 Ap 3 + ™ 2 — p 4m In (-^-=) - tt(R<) , (34) 

1714 \R C + J 

where f2(P<) is a constant for fixed (P, T, Ax). 

It is interesting to notice from Eq. fl3~2|) that if Ap 3 > p c the vortex line loses its 
metastable character. Indeed, when p 3 = p c one has R > = P<: the energy barrier disappears 
and the vortex core freely expands triggering phase separation. At P = 0, taking n — 1, 
Pim ~ 0.020A and the experimental values for the other variables, one gets P ~ 7.6A 
and He ~ 0.032 K which yields Ax cr ~ 1.3%. It is remarkable that this good agreement with 
experimental supersaturation values comes from an extremely simple model. The smallness 
of the obtained Ax cr fully justifies the use of the hollow core model to calculate Ax cr and 

To determine T* we proceed as in the preceeding subsection. Equation ( p2[ ) gives now 
the kinetic energy per unit length, and the expressions ( p3|) hold with r being the cylindrical 
radial coordinate. Again, the continuity equation in cylindrical coordinates can be formally 
integrated yielding 

1 r r 

U q (r,t) = -—r / S/j q (s,t)ds . (35) 



rp q (r,t) Jo 

The collective mass for this geometry is readily obtained: 



M(5) = 27cM(R c + 5) 2 \n[ 1 ^) , (36) 



where M. is defined in Eq. (|27|). In deriving the vibrational collective mass we have not 
taken into account any contribution from the circulation of the superfluid component, as only 
small radial displacements have to be considered if one supposes that phase separation is 
triggered by the expansion of the vortex core. Finally, we have for the crossover temperature 



k B T* = — 
2tt 



1 Ap 3 p 30 - R a/R 2 c _ 



M Rl In 

Taking ln^/Pc) ~ 1, one gets T* ~ 1 mK for Aa; ~ 1% at P = 
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IV. SUMMARY 



In this work we have calculated the crossover temperature and critical supersaturation of 
3 He- 4 He mixtures at very low temperatures using a functional-integral method in conjuction 
with either a density functional or a capillarity approximation description of nucleation 
configurations. We have found that the existence of vorticity in the mixture would explain 
the small Ax cr measured in experiments, but fails to reproduce T*, yielding values an order 
of magnitude smaller than experiment. Conversely, conventional nucleation seeds consisting 
in 3 He droplets do reproduce the experimental crossover temperatures but yield Ax cr values 
an order of magnitude larger than experiment. 

In spite of its failure, we feel that a detailed exposure of the merits and drawbacks of 
the present method might help find alternative roads to address this longstanding problem 
and identify any missing physical ingredient relevant to its solution. Besides, we think that 
our application of the FIA to this problem is especially simple and transparent, which adds 
a pedagogical value to it. 

A line of possible improvement might consists in introducing from the start a dynamics 
based on a zero rather than a first sound hydrodynamic description of the harmonic vibra- 
tions leading to T*. The density functional we have used reproduces the Landau parameters 
of the mixture that are relevant for such a zero sound description in the long wavelength 
limit; among other quantities, 3 He effective mass and the incompressibilities of both isotopes 
have been taken into account as fitting quantities to fix the density functional parameters. 
Moreover, the nucleation barrier is known to be well described by the new generation of static 



density functional calculations presented in Refs. £y.[M],[5l| It would be quite natural to intro- 



duce the zero sound dynamics much along the fruitful procedure followed in nuclear physics 
to describe nuclear collective excitations with density-dependent nucleon-nucleon efective 
interactions^. Such a program, which involves solving the Random-Phase-Approximation 
equations in imaginary time, has been formally developed to study a problem having some 
similarities to ours, namely spontaneous nuclear fission below the fission barriered. Here and 
there, a practical implementation of this scheme is quite a formidable task, especially in 
presence of vortices. 

Finally, we would like to point out that the experience gathered in the description of 
the nuclear giant monopole resonance indicates that for this precise mode, first and zero 



sound dynamics yield quite equivalent results, see chapter 13 of Ref. [44| and Ref. [46[ Thus, 
in the case of nucleation driven by 3 He droplets, as far as the problem of determining T* 
is amenable to study spherically symmetric, radial vibrations of the system, we expect that 
both dynamics yield rather similar T*. 
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APPENDIX: 



In this Appendix we give some hints about how we have derived the dynamical equations 
leading to T* in the case of 3 He- 4 He mixtures. We recall that by construction, the static 
density functional assumes that 3 He is in the normal phase, and all 4 He is superfluid. For 
this reason, the 'two fluid' hydrodynamics of superffuids, as described for instance in Ref. 
17, is of no applicability here. 

The first step is to write the static density functional of the kind proposed by Dalfovoli 



see also Refs |14],[42]) in a Galilean invariant forme 



£(p<i, T <iJq) = ^~ T 3 + ^~ T4 + 9s(pb, pi){ps r 3 - j 3 ) + gi(p 3 , pi){pi n - j 4 ) 

+ 9u(P3, P4)(ml p 3 r 4 + m\ p 4 r 3 - 2m 3 m 4 j 3 • j 4 ) + V(p 3 , p 4 ), (Al) 

where the kinetic energy r q and current j q densities are defined in Ref. ^9], and the functions 
g 3 , g 4 and g 34 are not identically zero if the atoms have an effective mass m* different 
from the bare mass (we have taken % = 1 thorough the Appendix). In the static case, the 
system is invariant under time-reversal and the currents vanishEil This allows to construct 
the g functions from the static density functional without having to reparametrize it to 
reproduce the basic thermodynamical properties of the mixture. Yet, casting the static 
density functional into a form like Eq. ( |A1[ ) leaves plenty of freedom, as it is easy to 
convince oneself. The grand potential density is then built from the energy density, the 
chemical potentials and particle densities: u = £ — J2 q P^qPq- 

One then applies to the 3 He and 4 He fluids a radial collective velocity u q (f,t) ('boost') 
deriving from velocity potential fields s q (r,t), u q (f,t) = Vs q (f,t) and takes into account 
how the particle, kinetic energy and current densities are affected by these boostsEa'EJ. If 
33 — 34 — in the static system, it is easy to check that the imaginary-time Lagrangian can 
be written as 



m qPqSq ~ ^Xlq Pq(^ S q f 



^G 34 (Vs 3 - Vs 4 ) 2 + [u(p q , T q ) - uj(p qm , T qm )\ , (A2) 



where \ G 34 = m\ m\ p 3 p 4 ^34. The imaginary-time Hamiltonian Tt is readily obtained from 

^ = ^m q Pq{^Sq) 2 + ^G U {VS 3 - VS 4 ) 2 - [u}{Pq,T q ) ~ Uj(p qm , T qm )\ , (A3) 

q 1 
and then the Hamilton equations 

5H 

m q s q = -— (A4) 

yield, respectively, the continuity and equations of motion: 

m 3 p 3 = _ m3 v(p 3 M 3 ) - V[G 34 (n 3 - u 4 )] 

m 4 p 4 = -m 4 V(p 4 M 4 ) + V[G 34 (m 3 - u 4 )] (A5) 
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m 3 



m 4 - 



du 3 
~dT 

~d7 



_ r i . 2 , 9g 34 



«3 - «4j - 



(tt 3 - M 4 ) 2 - — 



(A6) 



The above equations are quite similar to those obtained in Ref. [46] to describe nuclear 
collective vibrations, and at this point it is appropriate to recall the last paragraph of Sect. 
IV. Linearizing Eqs. ( |^5|JA6|) is straightforward and proceeds as outlined in Ref. 32. We 
do not write down the resulting equations, which are more involved than those given in 
the above reference because here we have kept in the equations the terms arising from 
the possible existence of a g 3 4 contribution. Rather, we finish the Appendix with a short 
discussion about possible choices for the g functions. 

The more obvious and simple one is to take #4 =1(734 = and g 3 p 3 = 1/(2777.3), where 
m^(p 3 , p 4 ) is the 3 He effective mass in the mixture@@. This leads to linearized equations 
to determine T* which are extremely cumbersome but still manageable!!. They have also 
been employed to obtain the density functional results presented in this work, and consistent 
with this choice is to write the collective kinetic energy in the capillarity approximation as 
indicated by Eq. (H). 

The more general expression g 3 p 3 + m\p±g 3 ± = 1/(2777,3) holds if one uses the accurate 
parametrizatiorSEl 



1 



2m* 3 



2m 3 



1 _ EL - £L 

PZc Pic 



(A7) 



where the fitting constants p 3c p^ c can be found in the above references. It can be checked 
that setting g± = and splitting the above relation into 



#3(P3,P4) 
#34(P3,P4) 



2777,3 
1 



El 
pic 



2777,37774 



2 | 2p 4 

P3c PZcPic 

P4 __2_ 

P4c 



pi 



(A8) 



the good fit to the experimental 777,3 ^ s retained, and 



7774 



7773 
7774 



P4_ 

pic 



+ 



2p 3 



P'icPic pAc 



P3 



(A9) 



which yields 7774/7774 = 1.31 for one 4 He impurity^ in normal liquid 3 He at saturation. This 
value is in agreement with the experimental ont 
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FIGURES 



FIG. 1. Schematic barrier Af2(<5) and inverted barrier well. 

FIG. 2. Schematic figure to illustrate the calculation of the collective mass. 

FIG. 3. T* (mK) as a function of Ax(%). Dashed lines are the results obtained within the 
capillarity approximation, and solid lines are the results of a density functional plus FIA calculation. 
The upper solid and dashed curves correspond to P = 1 bar, and the the lower ones to P = 0. 

FIG. 4. Density profile of the critical nucleation configuration obtained using either a density 
functional method (thick lines), or the capillarity approximation (thin lines). The values correspond 
to P = and Ax = 16% 
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